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Abstract 

We develop a gradient-flow theory for time-dependent functionals defined in ab¬ 
stract metric spaces. Global well-posedness and asymptotic behavior of solutions are 
provided. Conditions on functionals and metric spaces allow to consider the Wasser- 
stein space and apply the results for a large class of PDEs with time-dependent 

coefficients like confinement and interaction potentials and diffusion. Our results can be 
seen as an extension of those in Ambrosio-Gigli-Savare (2005)[2] to the case of time- 
dependent functionals. Eor that matter, we need to consider some residual terms, time- 
versions of concepts like A-convexity, time-differentiability of minimizers for Moreau- 
Yosida approximations, and a priori estimates with explicit time-dependence for De 
Giorgi interpolation. Here, functionals can be unbounded from below and satisfy a type 
of A-convexity that changes as the time evolves. 
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1 Introduction 


We consider the gradient flow equation 


u'it) = —V£{t,u{t)), t > 0, (1.1) 

m(0)=mo, (1-2) 


where £ : [0, cxo) x X —)■ (—cxo, cxo] is a time-dependent functional and (X, d) is a complete 
separable metric space. Our aim is to construct a general theory in metric spaces that can 
be applied for PDEs with time-dependent coefficients. In fact, with this theory in hand, we 
obtain global-in-time existence and asymptotic behavior of solutions in the Wasserstein space 
for a number of PDEs with density of internal energy U, confinement potential V and 
interaction potential W depending on the time-variable. That space consists of probability 
measures on Q with finite second moment endowed with the so-called Wasserstein metric 
d 2 (/U., u). Here we will focus on the whole space 

Gradient flows theory has been successfully developed for the case of time-independent 
functionals £{u) in general metric spaces (X, d) (see [7],[8],[2],[3],[15]). Two basic tools 
in the theory are the concept of curves of maximal slopes (see [8],[15]) and a time-discrete 
approximation scheme (see [7],[3]). The latter is based on the implicit variational scheme 

G Argmin{;^d^((7])“\r;) -hi^(n)}, (1.3) 

vgx 2 r 


where r > 0 is a time step. Notice that (1.3) consists in finding minimizers for interactive 
values of the Moreau-Yosida approximation S’t-{u) := + £{v)} of £ in 

(X, d). Speak generally, basic hypotheses assumed on £ are lower semicontinuity and some 
type of convexity and coercivity (see [2]). Eor the analysis of PDEs as a gradient flows, a 
suitable metric space is ^^2 in which the above theory has demonstrated to be particularly 
very fruitful. The idea of using the above discrete scheme in ,^^2 goes back to the work 
[11] for the linear Eokker-Plank equation and [20] for the porous medium equation. Sub¬ 
sequently, several authors extended this approach to a general class of continuity equations 
(see [2],[1],[6]) with velocity field given by the gradient of the variational derivative of a 
time-independent functional, namely 


du 

dt 


div 



in (0, -foo) X 


(1.4) 


where £ is the free energy associated to PDE dealt with. Under some basic assumptions, 
they considered £ with the form 


£[u] 


U(u(x)) dx -f 


u{x) V(x) dx + 


1 

2 


W{x — y) u{x) u{y) dx dy, (1.5) 


where U : R"*" ^ R is the density of internal energy, V : R'^ —)■ R is a confinement potential 
and W ; R'^ —)■ R is an interaction potential. The functional (1.5) has the classical form given 
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by the sum of the internal energy, potential energy and interaetion energy funetionals that is 
verified by a wide number of physieal models. Beside existenee of global-in-time flows, 
the literature eontains results on uniqueness, global eontraetion, regularity, and asymptotie 
stability of solutions (see e.g. [2]). We also quote the paper [4] where a ID non-local fluid 
mechanics model with velocity coupled via Hilbert transform was analyzed by using gradient 
flow theory in ^ 2 - 

In [13], the authors dealt with nonlinear diffusion equations in the form 

dtu - diw{A{V{f{u)) + uVV)) = 0, 


where A is a symmetric matrix-valued function of the spatial variables satisfying a uniform 
elliptic condition and /, V are functions satisfying suitable hypotheses. They also analyzed 
the contraction property for solutions. 

On the other hand, from a theoretical and applied point of view, it is natural to consider 
a time-dependence on the coefficients of some equations. For instance, a version of the 
stochastic Fokker-Plank equation (the one considered in [11]) is 


dXt = —'VV{t, Xt)dt + ^/2K{t)dBt, 


( 1 . 6 ) 


where the term is known as the diffusion coefficient and Bt stands for the classical 

Brownian motion. For (1.6), it is well-known that the law of processes is modeled by the 


PDF 


dtu = K(t)Au + V ■ (Vl^ (t, x)u). 


Another example is the version of the Mckean-Vlasov equation [24] 


dXt = b(t, fit, Xt)dt + \/2K{t)dBt, with b(t, fi, x) = —VW{t, •) * /i, 


where fit is the law of the processes Xt that obeys the PDF 

dtU = kAu — V ■ {b{t, u, x)u) 

with K depending on the time t. The term b{t, u, x)u corresponds to an interaction between 
particles with time-dependent potential. 

For a bounded convex domain C and 0 < T < cxo, Petrelli and Tudorascu [21] 
considered the non-homogeneous Fokker-Plank equations 

Ut-Xx- {uXx'ipit, x)) - Ax(P(t, u)) = g{t, X, m) in D x (0, T) (1.7) 

with Neumman boundary conditions and nonnegative uq E L°°(f2) such that / uodx = 1. 
They proved existence of nonnegative bounded weak solutions by constructing approxi¬ 
mate solutions via time-interpolants of minimizers arising from Wasserstein-type implicit 
schemes. Let us point out that, when P{t, z) = K{t)z, the conditions in [21] require that the 
viscosity k is bounded away from zero, while here we allow k to be arbitrarily near zero (see 
Theorem 6.9 in subsection 6.2). 
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In [23], Rossi, Mielke and Savare analyzed the doubly nonlinear evolution equation 

d'lpiu^t)) + duS{t,u{t)) 3 0 in B', a.e. te(0,T), (1.8) 

where B is a separable Banach space, 0 < T < oo, and m(0) = Uq. They proposed a 
formulation for (1.8) in a separable metric space (X, d) that extends the notion of curve of 
maximal slope for gradient flows in metric spaces. Existence of solutions is proved by means 
of a time-discrete approximation scheme in (2f, d) defined as 

(/" e Argmin{r?/;(d((7”“\u)/r) + S{tn,v)}, (1.9) 

v&X 

where r is a partition for [0,T] and r = |t| is the time step. Among others, the authors 
of [23] assumed that 8 satisfies the chain rule, is locally (in time) uniformly bounded from 
below, and differentiable in the f-variable with the derivative satisfying the condition 

\dt 8 {t^ u)\ < C{ 8 {t, u) + d{u*^ u) + 2Co), (ITO) 

for some m* G X, where 

Cq = — inf 8 (t, u). 

t€[0,T],u&X 

In fact, functionals in [23] are the sum of two time-dependent functionals 8 i and 82 where 
81 is bounded from below and Aq- convex (uniformly with respect to f G [0, T]), and 82 is 
a dominated concave perturbation of 81 . For a bounded domain Vt and uq G ifo(f2), using 
the above approach for dil}{u'{t)) = u'{t) (gradient flow case), they also analyzed (1.1) with 
du 8 {t,u) = —^u + F'{u) — l{t) in the L^(f2)-metric. Theseresults were improved in [18] by 
considering more general dissipation ip. Moreover, in [18] the condition (1.10) was relaxed 
to \dt 8 {t, u)\ < C 8 {t, u). We also refer the reader to [19, 22] for stability results for doubly 
nonlinear equations in Banach spaces. 

Since our functionals are not bounded from below neither satisfies a estimate like (1.10), 
we can not to apply the theory from [23] and [18]. Here we assume the conditions El, 
E2, E3, E4 and E5 given in Section 2 (see pages 6 and 7). Notice that E4 gives some 
local-in-time control from below for 8 but allows it to be unbounded from below at each 
t > 0. In (1.10) it is required some control of the time-derivative of 8 in terms of the 
functional itself. Instead of such estimate, we work with a condition on the difference of 8 
in two different times (see E3). In order to recover the contraction property, inspired by the 
convexity used in [2], we propose a type of A-convexity that changes as the time evolves (see 
E5). Thus, functionals could “lose convexity” in a such way that the approximation between 
two solutions for large times still holds, because the contraction property depends only on 
the mass accumulated by A, i.e. f* \ {s)ds. In general the function A(f) can be unbounded 
both from above and below in [0, cxo) but, for the contraction, it is assumed to be continuous. 
In Section 6, we show how to extend results for the case of 8 {t, u) having a more general 
density of internal energy U{t,u) and viscous term —A^{P{t,u)) (see Theorem 6.10 and 
Remark 6.11 in subsection 6.3). There, the conditions on potentials prevent £^(f, p) to satisfy 
E3. In the case P{t, z) = K{t)z, the diffusion coefficient k is non-increasing. This condition 


4 


is necessary in order to have the uniform limit of the approximate solutions (4.8) in all finite 
interval [0, T]. 

Another application of time-dependent gradient flows appears in the context of pursuit- 
evasion games. Jun [12] considered gradient flows in suitable playing fields and investigated 
existence and uniqueness of continuous pursuit curves that are downward gradient curves for 
the distance from a moving evader, i.e. a time-dependent gradient flows. In fact, his result 
works well in CAT(A')-spaces (with K = 0) that are complete metric spaces such that no 
triangle is fatter than the triangle with same edge lengths in the model space of constant 
curvature K. Also, he assumed that £{t, u) is Lipschitz in t, locally Lipschitz in u, and Aq- 
convex for alH > 0 where Aq is a fixed constant (i.e. Ao-convex uniformly in t). Another 
basic hypothesis used by him is that S’t^riu) given by 

S’t^riu) := mf{-^(f{u,v) + £(t,v)} 

vex ZT 

is Cr-Lipschitz in t, for all m e X = C'AT(O) and r > 0, where (7 > 0 is a constant. For 
the time-independent case £{u), we refer the reader to [16] for X = CAT{0) (see also [2]) 
and [14] for a geometric approach in X = CAT{K). 

In this paper we follow the program in the book [2] that contains a relatively complete 
gradient flows theory in general metric spaces and its applications for the non-vectorial space 
^2 by using optimal transport tools. So, our results can be seen as an extension of those in 
[2] in order to consider time-dependent functionals. For that matter, due to time-dependence 
of £, we need to handle some residual terms (see e.g. (4.14) and the estimate (4.19)) and to 
consider time-versions of concepts like A-convexity (see E5) and interpolation functions as 
(4.3) to (4.7). One of these functions is the interpolation (4.4) that corresponds to the time- 
dependent convexity parameter A(t). Thus, some adaptations from arguments in [2] made 
here is not a straightforward matter and involves certain care. Also, the time-differentiability 
of the minimizer for the Moreau-Yosida approximation of £ needs to be analyzed (see Propo¬ 
sition 3.4) and, in order to get the convergence of the approximate solutions, a priori esti¬ 
mates with explicit dependence on the t-variable are performed in Proposition 5.2 for which 
the aforementioned condition E3 plays a key role. 

The plan of this paper is the following. In Section 2 we recall some concepts such as 
proper functional and local slope, and some results on gradient flow theory in metric spaces. 
Also, we give the metric formulation for (1.1)-(1.2) and the basic assumptions for the func¬ 
tional £. In Section 3, we construct the approximate solutions, provide some properties for 
the minimizer of the Moreau-Yosida approximation, and give estimates for approximate so¬ 
lutions. In Section 4, we derive a priori estimates for the approximate solutions and show 
their locally uniform convergence in [0, cxo). In Section 5, we show that the curve, which is 
the limit of the approximate solutions, is in fact a solution of (1.1)-(1.2) in the sense of Sec¬ 
tion 2 and obtain the contraction property for solutions. Section 6 is devoted to applying the 
general theory in the Wasserstein space for PDFs with time-dependent functionals as those 
mentioned above. 
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2 Metric formulation and implicit scheme 


Let {X, d) be a complete separable metric space and consider the functional £\ X ^ 
(—cxo, +cxo] . Recall that 8 is said to be proper whether there is mq G X such that 8 (mq) < cxo, 
and its domain is defined by 


Dom(£^) = {m e X : 8{u) < cxo}. (2.1) 

Thus, a functional 8 is proper when Dom(£^) ^ 0. Let /+ and /” denote the positive and 
negative parts of an extended real-valued function /. The following concept is crucial in the 
theory of gradient flows, and we will use it for the case of time-dependent functionals. 

Definition 2.1. Let be a proper functional in a metric space X. The local slope \d8\ of 8 
at the point m G X is defined as 

|a£|(«) = limsup<L4zXX, (2.2) 

d{u,v) 

In what follows, we recall a technical lemma that will be useful in our calculations. 


Lemma 2.2 ([2, Lemma 2.2.1]). Let 8 : X ^ {oo, cxd] be a functional such that there is 
r* > 0 and u* E X with 


Sr*{u*) 


inf 

v&X 



d‘^{v, u*) \ 
2t* J 


> —oo. 


Then 

and 


> S‘r*{u*) - d^{u*, m), for all 0 < T < T* and u E X, 


— T 


4r*r / df{u^v) 


d {u,v) < — - 8{v) + 


T — T 


2t 


S’r*(u*) H- d‘^(u*,u)] . 

T* — T J 


In particular, the sub-levels of the map v ^ 8{v) + 


are bounded. 

Zt 


2.1 Metric formulation 


Let 8 : [0, -foo) x X ^ (—oo, +oo] be a time-dependent functional. It is well known 
that the problem ( 1.1 )-( 1 .2) admits a metric reformulation by using the concept of local slope 
(see [23]). This is given by the variational inequality 

2(f((,„(())) < S,£•(«,«(«)) - (i)P(«W) -(2.3) 


where \d8{t) \ stands for the local slope of the functional u -E- 8{t,u), for each fixed f > 0, 
and 


-u'l(f) lim 

S^t 


d{u{s),u{t)) 
|t — s| 
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stands for the metric derivative of an absolutely continuous curve u. 

Below we state the principal assumptions on the family of functionals £{t, •) on X, for 

t G [0, cxo): 

El.- For each f > 0, •) is proper and lower semicontinuous with respect to the metric 

E2.- The domain of the functionals, D ;= Dom(£^(f, •)), is time-independent. 

E3.- There exist u* e X and a function f3 : [0, oo) —)■ [0, oo) with f3 e oo)) such 

that, for each m G D, the function t ^ £{t, u) satisfies 

\8(t,u) — £{s,u)\ < [ (3{r) dr{l + (f{u,u*)). (2.4) 


Note that if the condition (2.4) is valid for some u* E X then it is in fact valid for all 
u* G X. Also, for each m G D, the function t ^ £{t, u) is differentiable a.e. in [0, oo) and 
its set of differentiability points may depend on u. 

Now we are ready to give the notion of solution for (1.1)-(1.2) that we deal with. 

Definition 2.3. Let uq E X and £ : [0, +oo) x X ^ (— oo, +oo] be a functional satisfying 
the assumptions El, E2 and E3. We say that an absolutely continuous curve u : [0, +oo) ^ 
X is a solution for (1.1)-(1.2), if m(0) = uq, the function t -E £{t,u{t)) is absolutely 
continuous, 

Iti'l, |8£(-)|W-))eiL(|0.«))). (2.5) 

and the variational inequality (2.3) holds true. 

2.2 Implicit variational scheme 

We start by recalling the Moreau-Yosida approximation of £. For r > 0 and f > 0, this 
approximation is defined as 


4,^(m) := inf{E(t,r,M;u)}, (2.6) 

v£X 

where the functional E(t, r, u] •) is given by 

Ei£T,u;v)-.= £i£v) + ^-^. (2.7) 

Zt 

Next, take a partition r = {0 = < • • • < < • • ■ } of [0, oo) with lim„_^oo = 

oo. Defining the step size Tn '■= one can construct the sequence 

f/; G Argmin{E(f!,^,r„,f/;-^r;)}, (2.8) 

v^X 


for a given family of initial data E X. 
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Since the eonvergenee results are locally in time, we can fix T > 0 arbitrary and analyze 
the convergence in [0, T]. In order to analyze rigorously the problem of minimization (2.8), 
we give two additional assumptions that will allow to obtain uniqueness and a nice behavior 
of the minimizers. 

E4.- For each T > 0, there exist a u* E X and r*(T) = r* > 0 sueh that the function 
t is bounded from below in [0, T]. 

E5.- There is a function A : [0, oo) —)■ R in oo)) such that: given points u, vq, vi G 

X, there exists a curve 7 : [0,1] —)■ X satisfying 7 ( 0 ) = vq, 7 ( 1 ) = vi and 

E(f, r, u] 7(s)) < (1 - s)E{t, r, u] Vq) + sE{t, r, m; Vi) - ^ - s)d‘^{vo.Vi), 

(2.9) 

for 0 < r < ^ and s G [0,1], where = max{— inftg[o,r] A(t), 0}. 

Remark 2.4. Note that by Lemma 2.2, for eaeh 0 < r < r* and u E X, we have that the 
function t —)■ St^T-{u) is bounded from below in [0, T]. In view of the assumptions E4 and E5 
we assume by technieal reasons that r* < min{-2—, 1 }. 

'^T+1 

Remark 2.5. In E5, we eonsider the existenee of eurves 7 : [0,1] —)■ X for all wq, vi E X and 
not only for elements in the domain D. This will be neeessary for the applications in Section 
6 where we will use the eoneept of generalized geodesies in the Wasserstein space 
These eurves exist independently of the funetionals that we will analyze in that seetion. 


3 Construction and properties of the implicit scheme 


In this seetion we provide some results about the sequence defined in (2.8). They can be 
seen as extensions of some results in [2] to the case of time-dependent functionals. We start 
with the following preliminary result. 


Lemma 3.1. Suppose El, E4 and E5 and letuEX,0<t<T, 0<t< -X, Then, the 
minimization problem 


min{E(f, r, u] u)} 


has a unique minimizer u\. 


Proof. Let Vn E D be a minimizing sequenee, i.e., lim E(f,r,M;f„) = St^r{u). Given 

n—^oo ’ 

m, n G N, by the eonvexity property E5, there is a curve 7 : [0,1] —)■ X, with 7 ( 0 ) = Vn, 
7 ( 1 ) = Vm, and 


< E(f,r,M; 7 ( 1 / 2 )) 

< iE(f, r, u; Vn) + ^E(f, r, u] Vm) 


T 


-1 


+ 


d ipni '^m) ■ 
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Thus 


- — ^ cf{vn,Vm) < (E(t,r,M;Un) - . 

It follows from the above estimative that Vn is a Cauchy sequence in X and hence it converges 
to some u\ E X. From the lower semicontinuity, we get that is a minimizer of the 
functional E(f, r, m; ■). The uniqueness follows from E5 and is left to the reader. □ 

In the next lemma we show that and u\. depend continuously on (r, t, u). 

Lemma 3.2. Assume the properties El to E5. Then, the following statements hold true: 

a) The map (r, t, u) G (0, r*(T)) x [0, T] x X —)■ ^^{u) G M is continuous. 

b) The map (r, t, u) G (0, t*{T)) x [0, T] x X ^ G X is continuous. 

Proof. We start with item a). Let {Tn,tn,Un) be a sequence converging to {ro,to,uo) in 
(0, t*{T)) X [0, T] X X. Denote by Vn = the minimizer of E(f„, r^, •) given in 

Lemma 3.1. It follows that 

limsup(ft„,r„(Mn) = lim sup E(f„, r„, 

n—^oo n—>-cx) 

< limsupE(t„,rn,Mn;w) 

n—>-CX) 

= E{to,To,uo;v), 


for all f G X. Taking the infimum in the right hand side, we obtain Ymisnp St^^rr,{un) < 

n—>-oo 

^*0,7-0 (wo)- In view of Lemma 2 . 2 , the sequence Vn is bounded. So, we can estimate 


liminf > Imvini \ ^{d{un,uf) - d{vn,uo)f + £{tn,Vn) 

n^oo n^oo I ZTn 


= lim inf 


d {Uni ^o) d{Un, uf)d{UQ, vf) Tq 2 


+ 


2 Tfi Tji 2 ,T()Tji 

+ E{to,To,Uo;Vn) + {£{tn,Vn) - (to,'^n))} 

> \imM\s’to,ro{uo)- [ l3{r) dr{l + d^{u*,Vn)) 

— <^to,To{tlo), 

and thus lim (^t^^TAun) = <^to,Toiuo)- For item b), note that 

n—>-CXD ’ ’ 

T./, . ^ ^ ^ _ f d‘^{Uo,Vn) d‘^{Vn,Un) 

'^nj ^tn,rn\^n) I 

V ^0 Tn 

-\-(^S (tf), Vn) S(tnj Vn)). 


d%Vn,Uo) 
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Using E3, the boundedness of Vn and the eonvergenee (r„, ^ (tq, to, uq), we get 

lim (E{to,To,Uo;Vn) - 4„,T„(Mn)) = 0, 

n^oo 

and, by the item a), 

lim E{to,To,Uo;Vn) = S‘t^,ro{Uo). 

n—>-cx) 

It follows that Vn is also a minimizing sequenee for S’to,To{uo) and then, by the same argu¬ 
ments in the proof of Lemma 3.1, it converges to (Mo)to’ required. □ 

Because of Lemma 3.1, for each family of initial data f/° G X associated to a partition 
r of [0, -|-cxo), we have that the sequence (2.8) is well-defined for each n G N such that 
< T -f r*, if r* < In what follows, we give some estimates for the minimizer of the 

Moreau-Yosida approximation (2.6). These will play an important role in the convergence 
of approximate solutions. 

Lemma 3.3. Assume that £ satisfies the properties El to E5. Let 0<r<r*, 0<t<T 
and M G D. Then 


,u*) — d^{u,u*) < e 




+ r 


r 


(3.1) 


for all e > 0. Ifr < r*/8, we have 

d^{v!'^'^,u*) < At* + 

+ Ad^{u,u*). 

Proof. We have that 

p(^,t+T 


ft+T* 


f3{r) dr{l -\- d^{u,u*)) — inf Sn^T-*{u*)) 


0<r<t+T* 


(3.2) 


dfiu^^fu*) - d\u,u*) = -2d{ul+fu*){d{u,u*) - d(<+",M*)) 

— {d{u,u*) — d{v!'^'^ ,u*)Y 
< 2d{u\+\u*)d{u\+\u) 


< e———— + T———^—- 


(3.3) 


< 2e{£{t + T,u) - £{t + T,u.f'^)) + T 

< 2e{£{t + T,u) - S’t+r,r*{u*)) 

^ d‘^{u^fi'^,u*) ^ df{v}fi'^,u*) 




(3.4) 


e 

Notice that we have already obtained (3.1) in (3.3). Now, choosing e = ^ in (3.4), we get 

]^df{u^^'^,u*) < T*^£{t,u) + j^ (3{r) dr{l +d‘^{u,u*)) - (g’t+r,T*{u*)) 


2t 


+ii^{u,u*) -\ M*), 




which implies (3.2) when r < r*/8. 


□ 
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The next result gives a time-differentiability property for 

Proposition 3.4. Assume El to E5. For 0 < r < ^, the function r St+T,T{u) is locally 
absolutely continuous in (0, ^hen is differentiable almost everywhere in that interval. 

For each m G D, assume further that the set of differentiability points oft^ S(t, u) does 
not depend on u (e.g., when t ^ S(t, u) is differentiable). Then 

^^t+rAu) = dtSf + T, <+") - ( 3 . 5 ) 

in the set of differentiability points. 

Proof. Let 0 < Tq < ri < ^. Reealling that minimizes + tq, tq, m; •), and using 
E3, we have that 


<^t+Ti ,ri A) — (Ft+ro ,toA) < (t + Ti, ) — 8 {t + Tq, ) 

(3-6) 

< [ P(T)dr(l + d‘‘(u‘+^\u’)) + '^ — 

Jt+To 

Similarly, but now using it follows that 

Ft+n ,Ti A) - Ft+ro ,ro A) > F{t + Ti, )-8{t + To, ) 

(3.7) 

ZTiTo 

> ^ — —d^A,uAA-[ Ar)dr{l + d^AlA'^^*))- 

‘^'Tl'To Jt+To 

Notiee that (3.2) allows us to estimate the terms dA*,uAA dA^u^A^ i = 0,1, by 
an expression independent of r, which gives the absolute continuity in each compact in¬ 
terval of (0, ^]. Now take a point r G (0, where the derivative of r —)■ ^i+r,r exists. 
Considering lateral limits, the equality (3.5) follows by using estimates (3.6)-(3.7) and that 
dAA'^^ 'A!fA —)■ 0 as Tfc —)■ r (see Lemma 3.2 b)). □ 

As a consequence, we have the following corollary. 


Corollary 3.5. Assume the same hypotheses of Proposition 3.4. Then, for m G D, we have 
the identity 


dfA-^Af'^) r dfA-:AA'^) 


2t 


dr = I dt8{t + r,uA'')dr + 8{t,u)—8(t + T,uA^)- (3.8) 


£+T\ 


2 r 2 


Proof. By integrating (3.5) from tq to r < ^, it follows that 




2 r 2 


dr = I dfSit + r, dr. 
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In view of the definitions of St,T{u) and u\., and since the above integrals are finite as Tq —)■ 0, 
the remainder of the proof is to show that St+TQ,TQ{u) — S{t, u) as Tq -> 0, for each fixed 
f > 0. In fact, note that 

^+ro,To(^) 


and so 


limSUp4+ro,ro(M) < £{t,u). 

ro^O+ 

Also, we can conclude from (3.9) and Lemma 2.2 that d{u, 
the lower semicontinuity of we get 


(3.9) 

—)■ 0, as To —)■ 0. Using 


£{t,u) > 

lim sup 

^t+TO: 

TO 


TO—>-0 



> 

lim inf I 


“to 


TO—>-0 

> 

S{t,u), 





/9(r) dr{l + d'^{u*, 



as desired. 


□ 


Remark 3.6. In the last proof, we have showed in particular that —)■ m as r —)■ 0, when 

M G D. 

Now, we recall a discrete Gronwall lemma. 

Lemma 3.7 ([2] Lemma 3.2.4). Let A, a G [0, oo) and, for n>l, let an,l3n £ [0, C)o) satisfy 


n 


ttn < A + a fdjttj, Wn > 1, with m 
i=i 

Then, denoting B := A /(1 — m), 6 := a /(1 — m) and {3o 


SUpQ!/9„ < 1. 
nGN 

: 0, we have that 


On < n > 1 . 


The variational scheme (2.8) will be the base for constructing approximate solutions for 
(L1)-(L2). The below lemma can be seen as a version of [2, Lemma 3.2.2] for the case of 
time-dependent functionals and gives a first set of estimates in order to control approxima¬ 
tions. 

Lemma 3.8. Assume El to E5. Let t = {0 = < t]. <■■■< 1^. <■■■} be a partition of 

[0, oo), Tj := and |t| = sup ■ |r, |. ForT > 0 and t* < choose G N such 

J ^T+1 

that T G t^)- Suppose that there is a constant S' > 0 satisfying 

S(0, U°) < S and d\uf U°) < S. (3.10) 

Then, there exists a constant C = C{S,T,t*,S) >0 such that 

d\u*, Uf) <67, ^ ^ Ui-^) - Ui)) < 67, (3.11) 

j=i 6 ?' 

for all 1 < n < N and \ r | sufficiently small. 
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Proof. By the minimizer property of Ul and E3, we get 

j=i j=i 

+ J2 I dr{l + d\u*,Ui-^)), 


(3.12) 


for 1 < n < iV. Using the first estimate in Lemma 3.3 with u = Ul = Ui and 

e = Y’ obtain 




n 

Y^-d‘{u\Ui)--d‘{u\U=, 


- 1 ^ 


i=i 


< - £(0,(/°)- inf + 


0<t<T+T* 


d\u*,U^^ 




d^Ufu*) T 


* rt-T 


J=l 




+ 


/3(r) dr ) d‘^{u*, 


2 /+j -1 


y 


r-T+r* 


/9(r) dr. 

^ Jo 

Rearranging terms, it follows that 

d^u-, U") < 2 t-{s- _ («•)) + 2 (l + r- jf' /J(r) dr) S 

^{r) dr+ A^ + ^ I3{r) dr\d^{u\Uf) 


0<t<T+T* 
r-T+r* 


+2r* 


'0 


i=i 


< r4(^,T,r*,^)+4 5^/3,d2(«*,f/^), 

i=i 


(3.13) 


j+i 


for some constant A = A(S, T, t*, 8) > 0, where fdj '■= ^ + \ (dip) dr. By using an 
argument of absolute continuity, we have that max 4/9, < 1, for |t| small enough. Then, 

l<n<Af 

the first estimate in (3.11) follows by using Lemma 3.7 in (3.13). For the second one, we use 
(3.12) and observe that 

^2/ * TTn\ 


which is bounded. This concludes the proof. 


□ 
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4 A priori estimates 

It is well known that, under eonvexity hypotheses, the problem (1.1)-(1.2) admits a for¬ 
mulation based in a differential inequality. In faet, in the ease when X is a Euelidean space 
and the functional •) is A(f)—convex, the curve solution u{t) satisfies 

-vf + - vf + £{t,u{t)) < (4.1) 

for all V E X. Assuming the hypothesis of convexity E5, one can derive a discrete version 
of (4.1). In fact, for each fixed f > 0, we have (see [2, Theorem 4.1.2]) 

- ^d^{u,v) + \x{t)dF{u\,v) < S{t,v) - St,r{u). (4.2) 

It It I 

Now we define a set of interpolating functions that will be useful in the convergence of 
approximate solutions. In comparison with [2], the time-dependence of E generates new 
residual terms in the estimates and leads us to define the interpolations Tt and X-rit) in (4.3)- 
(4.4) below. The function X-r^t) is necessary in order to deal with the time-dependence on 
the parameter A. 

Let T = {0 = < • • • < < • ■ ■ } be a partition of [0, oo) and := 

Consider T > 0, X G N such that T G t^], and the following functions defined on 


the interval [0, T]: 

TAt) ■.= tA for te{tA\K], (4.3) 

XAt):=X{tA, (4.4) 

f — 

lAt) ■= -for t G (4.5) 

'^n 

dl{t-V) ;= {l-lAt))d\Ur\V) + lAt)d\U^,V), for t G (4-6) 

£At) ■■= ii-iAt))£{tA\K-^) + iAt)£iK,u^A, for t g (4.7) 

UAt) := K-\ UAt) ■■= UA for t G (4.8) 

Also, we consider (/^(O) = U.A^) ■— Th® functions in (4.8) are called approximate 
solutions for (1.1) corresponding to the data (7°. 

Taking A = UA “ = and v = V,we can rewrite (4.2) as 

~4(i; V) + V) + f,(() - £(TAt), V) < (4.9) 

Arit) + (1 - UT') - m, c'r')). 

for t G A], where 

A-it) := (1 - -£K,u;)) - (4.io) 

With this notation, we have the next estimate. 
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Lemma 4.1. Assume El to E5. For a partition r with |t| < t*, define the residual term 


:= {l-lr{t))d{U^{t),U,{t)). (4.11) 

We have that 

i44((; V) + V) - (AJ(«)<i(F,((), £/,(«)) +a;(()S?.(«)) d^(t-, V) 

+ f.(i) - £iTAt),V) < + 4^®? + (1 - 4)(£(«r'.e^) - £K,UA), (4.12) 

for all V G D and almost every point t G [0, T], 

Proof. A detailed proof for the case AT-(t) < 0 can be found in [2, pg.88]. Let us explicit the 
proof for AT-(f) > 0. For that, we can suppose that 

d{u:^,v)<d{u:^-\v), 

and estimate 


d^{U^{t), V) - dl{t- V) = {1- W) {d\U^{t), y) - d\U^{t), V)) . 

Thus, we have 


d\U^{t),V)-dl{t-V) 


> -d{Ul^,U^-^){d{Ul^,V) + d{U^-\V)) 

+l^d{Uf, U'^-^){d{U^-\ V) - d{U^, A")) 

+l^d{Uf,V) + d{Uf,V)) 

> -2d(t/;, Uf-^) ((1 - f")) 


which together with (4.9) gives the desired result. □ 

The next result is a slightly modified version of the Gronwall Lemma in [2, Lemma 4.1.8]. 
The proof is the same and we omit it. 

Lemma 4.2. Let x : [0, cx)) —)■ M a locally absolutely continuous function and let a,b,X & 
oo)) be such that 

+ 2\{t)x^{t) < a{t) + 2b{t)x{t) a.e. t >0. (4.13) 

(ajL 

For T > 0, we have that 


where a{t) = A(s) ds. 




ft 


+ 


x2(0) + sup / e2"(*)a(s) ds j +2 / e"^‘^|&(f)| cZf, 

te[o,T] Jo j Jo 
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4.1 More two interpolation terms 


In this subsection we consider two interpolation functions that depend on two partitions 
T and 77 of [0, oo) with |t|, [r/l < t*. So far, we have define two residual terms and 
in (4.10) and (4.11), respectively. Another one that we will work with is 


:= 2(1 - 4(()) |£(r,((),£,(()) - £(r,((), £/,(())] 

+ 2l,{t)[e{%{t),UAt))-e{TAt),Ur{t))], forie|0,r]. (4.14) 
Define also the interpolation function s) as 

s) = (1 - lrj{s))dl{t,U^{s)) + lr,{s)dl{t,Urj{s)). 

Taking in (4.12) a convex combination, with coefficients — and lri{t) for V = ILr^it) 
and V = Urj{t) respectively, we arrive at 

+ {K + \)d%{t,t) < 2 [A+d(I7,(f),[/,(f)) + A+d(I7^(f),f/^(f)) 

+ A" (f) {t) + A" (f) % {t) d-rr^{t,t)+ {t) 

+ ^r){t) + A.,. (t)^‘^(t) + (t)3)‘^(t) + G-rrjit) 

+ Gr^T-it). 

Now we can use Lemma 4.2 in the last inequality in order to estimate 

e“-Wd^^(f,f)< (d2(f/0,f/0)+ re2a.,(.) ^ (m+{s) + \-,{s)&1{s)) ds 

rt X 1/2 

+ j (G^t.(s) + G+^(s)) dsj 

+ (A+(s)d(f7^(s),f/,(s)) + A+(s)d(I7^(s),f/^(s)) 


(4.15) 


+ \{s)%{s)j ds, 

for all t > 0, where a-rriit) '■= Jq At-(s) + A^(s) ds. 


4.2 Convergence of the approximate solutions 

In this section, we deal with the convergence of the approximate solutions U-r and IJ_^. 
Using the minimizer property of Ul and direct calculations, one can obtain 




(4.17) 


/ /•* \ ^ ” 


e“-<‘)AJWd((7,W,C/,W)& <C(r)|rp5^(£(«i,C/r‘)-£«,l'i)), 

3 / j=l 

(4.18) 

for T e (t^-S t^], 0 < t < T, and n < iV. 

The above estimates give some eontrol on the residual terms and Next, we provide 
an explieit estimate for the residual term (4.14). This eould be useful to obtain eonvergenee 
rates of approximate solutions to the gradient flow solutions. Reeall the standard notations 

a Ab = min{a, b} and aV b = max{a, b}. 

Proposition 4.3. Assume El to E5 and the boundedness condition (3.10). Let r, 77 be two 
partitions of [Q, + 00 ) with |r|, |r7| small enough as in Lemma 3.8. ForT > 0, choose N^K G 
N, such that T E ] H t^]. There is a constant C = C{T, S, r*, £) > 0 such 

that 


G+ (f) dt < C{\t\ + |r/| 


(4.19) 


Proof. Denote f”]. For t\, let ki be the greatest integer satisfying < t\. If 

= t\. define Ji = J.^. Otherwise, ehoose rii < iV as the greatest integer with the property 
and define Ji = U ■ ■ ■ U In both eases, we have Tr^iFf) = t’f^ and then 


bJl 


(1 - 4(f)) 


f‘Tr{t)\/Tr,(t) 


'Tr{t)ATri(t) 


(3(s) ds dt < / (1 — /x(^)) dt 


.TAtF)yTr,{AP) 


'Jl 


(3(s) ds 


< (kl + l^l)/ (3{s)ds. 

Jo 

If = t^, define Otherwise, take the greatest integer k 2 E N sueh that 

t^~^ < 111 the case t^ = define Otherwise, take the greatest 

integer 712 < N such that < t^ and define = 1”^+^ U ■ • • U Noting that 7^(f”^) = 
t’^ and Trit) > Ff^^, we get 


'Jl 


(l-4(f)) 


l'Tr{t)VTri(t) 


' Tr(t)ATr,{t) 


(3{s) ds dt < / (1 — 4(f)) dt 


.rrf?)yTn(fF) 


02 


< 


|i7l) 


f Ljj 


l-q 


4i+4r.,(0 

/)(s) ds. 


f3{s) ds 


Proceeding inductively, and adding estimates obtained in the process, we arrive at 


(1 - 4(f)) 


I'TriOVTrjit) 


' Tr{t)A%,{t) 


(3{s) ds dt < (|t| + | 77 |) 


fT+r* 


(3(s) ds. 


(4.20) 
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Analogously, 



I At) 


f‘TT{t)yTr,{t) 


'TT{t)ATr,{t) 


/9(s) ds dt < {\t\ + \r}\) 


r-T+r* 


(3{s) ds. 


Adding (4.20) and (4.21), we get 


T nTr{t)yTn{t) 


As) ds dt < 2(|r| + I77I) 


T+t* 


/3(s) ds. 


(4.21) 


(4.22) 


Jo J TT(t)ATr)(t) Jo 

Now, reealling (4. 1 1) and the property E5, and using the first estimate in (3. 1 1), for f G [0, T] 
it follows that 


GtAt) < 2{1-lAt)) As)ds{l + dAu*,UA{t)) 

J TT(t)/\Tr,{t) 

r'TT{t)yTrt(t) 

+2lAt) / As)ds{l + d^{u\llA{t)) 

JTr{t)/\%){t) 

rTT(t)\/Tr){t) 

< C{S,T,T*,S) / AAds. (4.23) 

JTr{t)ATri{t) 

Finally, we eonelude by integrating (4.23) over [0, T] and using (4.22). □ 

In the present seetion and in Section 3, we have obtained some properties and estimates 
for £{t,u) and the implicit variational scheme (2.8) associated to the problem (1.1)-(1.2). 
After doing that, we are in position for proceeding as in [2, pag. 91-92] and showing that the 
approximate solutions (4.8) converge uniformly in [0, T] as |t| 0. 

Theorem 4.4. Assume El to E5 and the condition 

lim d{UAuo) = 0, sup £^(0,(7°) = S < 00 , (4.24) 


for some uq G D. Then, the approximate solutions Ut and converge locally uniformly 
to a function u : [0, 00 ) —)■ X satisfying m(0) = uq. Moreover, u is independent of the family 
U^. 

Remark 4.5. In fact, the convergence of the approximate solutions is valid for uq G D. 


Proof of Theorem 4.4. The proof follows essentially the same arguments in [2] by taking 
care of the time-dependence. We give some steps for the reader convenience. By taking a 
suitable convex combination, we arrive at 

d^fJAt), UAt)) < ^dlAt^ t) + 3C(|t| -f |r7|). 
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So, using (4.16), (4.17), (4.18) joint with Lemma 3.8, and the estimate (4.15), we obtain 

< (d\U«,U^) + C{\T\ + \r,\)+ ' 

+C(|t| + |r;|). 

We eonelude the eonvergenee by using Proposition 4.3 and the completeness of the space 

X. □ 

5 Regularity 

In this section we show that the function u obtained in Theorem 4.4 is, in fact, a solution 
for (1 . 1)-(1.2) in the sense of Definition 2.3. For that matter, we need to show some regularity 
properties for u. We begin by recalling the De Giorgi interpolation. 

Definition 5.1. Let {U^)n be a solution for the variational scheme (2.8), defined for t” < 
T + T*. Define the De Giorgi interpolation 

UT{t) = + 5), for t e (fT~^ K] 6 = t — 

as the unique minimizer of the functional n G X —+ 5, 5, 

We have that the De Giorgi interpolation also converges locally uniformly to the same 
function u in Theorem 4.4. 

Proposition 5.2. Assume the same hypotheses of Theorem 4.4. There is a constant G > 0 
independent ofr such that 

d\UAt)Mt))<\T\c\l + j-^^^ /3(r)drj, (5.1) 

for all t G Thus, Ut converges to the function u given in Theorem 4.4 a.e. in 

[0, T]. Moreover, the convergence is uniform provided that the function (3 in E3 belongs to 

^te([0,oo)). 

Proof. Let X G N be such that T G t^]- First, we will show that the discrete solution 

satisfies the inequality 

£(C,f/;)<£(0,f/°)+G, (5.2) 

for some constant C independent of r. In fact, by using E3 and the minimizer property (2.8) 
of Uf, we obtain 

Uf) < Uf-^) + r /3(r) dr{l + d{uf 
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Recall that d{u*,U^ is bounded by a constant C that depends on T and is independent of 
T. Proceeding inductively, it follows that 


<£{0,U^) + C 


fT+r* 


/3(r) dr, 


from where we get (5.2). Now, estimate (3.2) in Lemma 3.3 and (5.2) give 

d{UAt),u*) < 4T*{£{t:^-\U^{t))+C{l + d\u*,U^m) 

< 4r*(£(0,(7°) + C'), 

for t e f"]. Taking S = t — for t G t^], we can estimate 




26 


< 


26 ^ 

dHUrit),UAt)) 

2t„, 


+ £{tAUAt)) 


+s{t,uAt))-SK,uAt))) 


Rearranging terms and using E3, it follows that 


26 2t„ 


7r,-7^] UAt)) < Ut - d\U,(t), C/,(i)) + / m dr 


26 2t, 


n 

X [2+ d‘^{u\U^{t)) +d:\u\Urit))] . 


Recalling that £{0, f/°) < S and using (5.2), we obtain (5.1) and then the convergence of 
Urit) to u{t) in the set of Lebesgue points of /3. □ 

Before proceeding, let us recall a well-known estimate for the slope \d£{t)\. Recall that 
stands for the minimizer of E(f -|- r, r, u; •). Then G Dom{\d£{t + r)|) and 

\dS{t + r)|«+") < (5.3) 

T 

Under the convexity hypothesis E5, we have that the local slope \d£{t) \ is lower semicon- 
tinuous and 

\d£{t)\{u) = sup + \At)d{u,v)\ . (5.4) 

v^u V d{u,v) 2 J 

The next lemma will be useful to show -regularity for functions with a certain type 
of control in their variations. 
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Lemma 5.3. Let T > 0 and f,g,l3 G L^([0, T]) be such that 


\f{t) - f{s)\ < {g{t) + ^(s))|t - s| + 


(3{r) dr, 


for s < t. Then f G W^'^{[h,T — h]), for all 0 < h < T/2. 


Proof Since the funetion t —)■ /3(r) dr belongs to T]), we have the difference 

quotient property 


sup 


pT-h 

2 


— fiir) dr 

Ih 

h Jt 


dt < oo. 


(5.5) 


Using the notation 


we obtain 




fit + h)- fit) 
h 


rT-h 


rT-h 


\^hif)it)\dt < 


git)+git + h) 


ft+h 


< 2 


1,1 


+ 


pT-h 

2 pt-\-h 


— / fiir) dr 

Jh 

h Jt 


/3(r) dr 
dt, 


dt 


whieh gives the desired regularity by employing a differenee quotient argument. □ 

Now we are ready to show that the limit u in Theorem 4.4 is a time-dependent gradient 
flow in the sense of Definition 2.3. 


Theorem 5.4. Assume El to E5. The limit u : [0, oo) —?■ X in Theorem 4.4 is locally 
absolutely continuous and its metric derivative \u'\ belongs to oo)). Moreover, if the 

function t ^ Sit, u) is differentiable for u G D, its time-derivative is upper semicontinuous 
in the u-variable (with respect to the metric), and the property 


, I , 7/ ^ V n ^ r • f ^i^ni tin) ^i^i tin) ^ o cfj. \ 

tn i t, diUni u) ^ u as n ^ OO ^ lim ml- > OtS (t, u) 


tn t 


(5.6) 


holds true, then the function t ^ Sit, uit)) is absolutely continuous and satisfies the identity 


Sit,uif)) — Sit),uit))) = / dtSis,uis)) ds -/ |m p(s) ds-/ |c)i^(s)f (m(s)) ds. 

Jo 2 do 2 do 

(5.7) 

In particular, u is a solution for ( 1.1 )-( 1.2). 


Remark 5.5. Definition 2.3 does not eontain (5.6). Note also that this assumption is used to 
prove (5.7) and, in faet, is not neeessary to obtain the absolute continuity of t — Sit, uit)). 
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Proof of Theorem 5.4. Let T > 0 and denote by 


\Km 


Tn 


(5.8) 


the discrete derivative of UT{t) in each interval ^t]- Lemma 3.8, we have that 

\K\‘^{s) ds < C, (5.9) 

for each t G [0,T]. Thus, we can extract a sequence such that \Tk\ —)■ 0 and \U!^J 
converges weakly in L^([0,T]) for some function m. Fix 0 < s < f < T and choose 
P = p(s) and n = n{s) G N with s G and t G It follows from (5.8) 

and triangular inequality that 

d{U^J.s),U^M) < l‘‘\K,\(r) dr. 

Letting k —j- +cxo, and using the weak convergence, we conclude that u is absolutely con¬ 
tinuous and \u'\ < m. Also, after a change of variables, we can employ the identity (3.8) to 
obtain 



1 

2 



|(/^P(r) dr + 



d^{UAr),UAr)) 

2r2 


dr 



dtS{r, Ur{r)) dr 


+8f),UA-S{r^,UA- (5.10) 


For the above subsequence, we have 


S(t,u(t)) < liminf£^(f, UtA^)) 

k^oo 


liminf£(7;,(f),f/rfc(f)), 

k^oo 


and so, using (5.3) and (5.4), we arrive at 


^ J dr + ^ j |c)£^(r)p(M(r)) dr + £^(f, M(f)) 

< r‘'V; 

- 1,^=0 V2 I ^ 2r2 

r'r^i.d) 


< limsup / 

fc—>-oo Jo 

< [ dtS{r,u{r)) dr + S{0,Uo), 
Jo 


dtS{r, UtAt)) dr + £:(0 ,m(0 )) 


(5.11) 


where, by convenience, we have chosen (7° = uq (recall that u does not depend on (7° —)■ uq). 

Notice that in particular sup 8 {t,u{t)) < oo. On the other hand, in view of [2, Lemma 

te[o,T] 
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1.1.4 a)], there exist an inereasing absolutely eontinuous funetion s : [0, T] —)■ [0, L], whose 
inverse t is Lipsehitz, and a eurve u : [0, L] ^ X sueh that |m'|(s) < 1 and u{t) = u{s(t)). 
Considering the funetion cp{s) = S{t{s),u{s)) and using (5.4), it follows that 

V9(si) - (^(S 2 ) < (|c)£^(f(si))|(M(si)) + A:^C) |S 2 - Si| 

+(1 + C^) f I3{t{s))t'{s) ds, 

J SI 

for Si < S2, where C = sup d{u*, u{s)). Replaeing the roles of si and S2, we obtain 

sG[0,L] 


|9?(si) - <yC(s2)| < {\dS{t{si))\{u{si)) + \dS{t{s 2 ))\{u{s 2 )) + 2Xj.C) |S2 - Si 

rs2 

+ (1 + C^) / {3{t{s))t'{s) ds. 

J SI 


By using Lemma 5.3, we ean eonelude that ip is absolutely eontinuous and then £{t,u(t)) 
also does so. It follows that £{t,u{t)) is derivable at almost every point t G [0,T]. Let 
to G [0,T] be a differentiability point of £{t,u{t)) for whieh the metrie derivative \u'\{to) 
exists. Taking tn i to, we get 


d 

dt 


{S{t,u{t))) > liminf 


£{tn,u{tn)) -£{to,u{tn)) 


lim inf 

n—>-oo 


tn to 

£{to,u{tn)) -£{to,u{to)) d{u{tn),u{to)) 


d{u{tn),u{to)) tn-to 

= dtS{to,u{to)) - \dS{to)\{u{to))\u'\{to). 


Integrating the above inequality, and using (5.11), we obtain (5.7). 


□ 


Corollary 5.6. Under the hypotheses of Theorem 4.4. There exists a subsequence of parti¬ 
tions Tk such that Erkif) defined in (4.7) converges to t ^ £{t,u{t)) in oo)), and 

therefore a.e. in [0, oc) (up to a subsequence), where u is as in Theorem 4.4. 

Proof. We only need to show that, for T > 0, the functions and Qt defined as 


Ut)-.= £{tl,U^), for fG(fr\C], 
and 

g^{t)-.= £(tl-\U:^-^), for te(tl-\t% 

converge to f —)■ £{t,u{t)) in L^([0,T]), as |t| —> 0. First, note that for each partition 
{0 = to < < • • • < = ^} of [0, T], we can bound the variation of as 

L N 

(5.12) 

1=1 n=l 


i-T+t* 


+c 


fi(s) ds, 


(5.13) 
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where C > 0 is a eonstant independent of r. By Lemma 3.8, the summation in the right 
hand side of (5.13) is bounded, and therefore the total variation of f-r in [0, T] is uniformly 
bounded. Analogously, the total variation of g-r in [0,T] is uniformly bounded. It follows 
from [10, Chap. 5, Theorem 4] that there exist a subsequenee Tk and funetions A, B G 
L^([0,T]) such that —)■ A and gT,^ B in L^([0,T]), as A; —)■ oo. Also, it is not hard 

to show that A = B > 8{t, u{t)) a.e in [0, T]. Now, the same argument used in the proof 
of (5.7) can be used in order to show the equality A = £{t, u{t)) a.e. t E [0, T]. In fact, if 
-E- Ain Li([0,T]),then 


^ J |M'|^(r)dr + ^ J \dS{r)\‘^{u{r)) dr + A{t) 

< liminfl i r-“V;.P(r)<ir+ A*"’* 


fc—)-oo 


2r2 


+£(r.,(f),L.J) 


< 


dtS{r,u{r)) dr + S{0,Uo), 


and, using (5.7), we are done. 


□ 


Remark 5.7. As a consequence, we have that the solution u{t) E Dom(|c)£^(f)|) for almost 
every point t E (0, oo). 


5.1 Contraction property 


Consider the condition 


E6.- The function \{t) is continuous. 


Having at hand the estimates obtained in previous sections, the contraction property holds 
if we assume E6. Here we only sketch its proof for the reader convenience. 

For \{t) continuous, the interpolation A,- defined in (4.4) converges uniformly to A, as 
|r| 0, in each bounded interval of [0, cxo). Recall the following technical lemma [25, 

Lemma 23.28]. 


Lemma 5.8. Let F = F(t, s) be a function [0, oo) x [0, oo) —?■ M locally absolutely contin¬ 
uous in the variable t and uniformly continuous in s, and locally absolutely continuous in s 
and uniformly in t; that is, there exists a nonnegative m E oo)) such that 

|F(f, s) — F(f', s)| < f m{r) dr and \F{t, s) — F{t, s')\ < f m{r) dr, 

Jt' J s' 

where m does not depend on s in the first inequality and on t in the second one. Then, 
the function 6{t) := F{t,t) is locally absolutely continuous and, for almost every point 
to E [0, oo), we have 


d 

dt 


t=tod{t) < limsup 

tpo 


F{to,f) - 5{to) 
t-to 


+ lim sup 


F{t,to) - S{to) 


t-to 


(5.14) 
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Integrating (4.9) from s to t with 0 < s < t < T and taking the subsequenee given in 

Corollary 5.6, we can pass the limit and use E6 in order to obtain the inequality 

-(f{u{t),V) — -d'^{u{s),V)+ [ + £{r,u{r)) dr < f S{r,V)dr. 

J s J s 

(5.15) 

Let u, V be two solutions given by Theorem 4.4 with initial data mq, Vq E D, respectively. 
Recall that, by Lemma 3.8, both curves u and v are locally bounded. Also, 

{d"^{u{t), v{s)) — d‘^{u{t'),v{s))\ < d{u{t),u{t')){d{u{t),v{s)) + d{u{t'),v{s))) 

<C{T) [ \u'\{r)dr. 

Jt' 

where C{T) = 2 sup {d{u{t),u*) + d{v{t),u*)). Similarly, one can show the local ab- 

0<t<T 

solute continuity in the variable s for the function F(t, s) = d‘^{u{t),v{s)). It follows that 
d‘^{u{t),v{s)) verifies the hypotheses in Lemma 5.8. Next, using (5.15), a direct computation 
gives 

—d‘^{u{t),v{t)) + 2X{t)d'^{u{t),v{t)) <0, (5.16) 

(ajL 

for almost every point t E [0, cxo), which implies 

d{u{t),v{t)) < ‘^"d{uo,Vo). (5.17) 

Remark 5.9. The time-dependent functional £ can be “weakly” convex (A(f) < 0) at a 
certain t = to. In fact, we could have A(s) ds < 0 and solutions distance themselves. 
However, according to the behavior of A(f), the convexity could be improved (A(f) > 0 
and Jg A(s) ds > 0) as t increases. In this case, we would recover the time-exponential 
approximation between the solutions u and v. 

6 Applications for PDEs in the Wasserstein space 

In this section we apply the theory developed in previous ones for time-dependent func¬ 
tionals associated to PDEs in the Wasserstein space. This space has a very nice geometric 
structure and is suitable to address gradient flow equations. 

We start by recalling some definitions and properties of that space. We denote by 
the set of Borel probability measures in with finite second order moment, i.e. 
/i E if /i is a positive Borel measure, 

yu(R'^) = 1 and M 2 {fi) := / |xp d^{x) < oo. 

We can endow 0 ^2{W^) with the weak-topology or the so-called narrow topology by consid¬ 
ering the following notion of convergence: 

/Xfc ^/i as/c ^ cxo lim / f[x)d^k{x)= / f{x)d^{x), (6.1) 

jRd J^d 
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for all / G where stands for the set of bounded eontinuous funetions. On 

the other hand, endowed with the Wasserstein distanee is a eomplete metrie space. 

This metric is defined by means of the Monge-Kantorovich problem and reads as 


d 2 (/i, = min i / \x - yf d'y{x,y) :-f e , (6.2) 

L./RdxIR'^ J 

where r(/i, z/) = {7 G x M'^) : 7(^4 x M'^) = 7(M'^ x B ) = In fact, 

there exists at least one probability measure in that reaches the minimum in 

(6.2). This is called the optimal transport plane and is supported in the graphic of the sub- 
differential of a convex lower semicontinuous function (see [26]). We denote by ^ 2 ,ac(R^) 
the set of probability measures in that are absolutely continuous with respect to the 

Lebesgue measure. If /i does not give mass to sets with Hausdorff-dimension less than d — 1 
(e.g. if /i G ■^ 2 ,ac(R^)), there exists a map tj) : —)■ that coincides with the gradient 

of a convex lower semicontinuous function, such that u = and the optimal transport 

plane 70 is given by the push-forward of /i via the map Idx i-e- 7 o = {Id X t;))#/i. Thus, 
we have that (see [26, theorem 2.12]) 

dl{f^,^)=[ \x-t''^{x)\‘^ dfi{x). (6.3) 

jRd. 

We also recall the concept of generalized geodesics [2]. 

Definition 6 . 1 . Let a, /xo, /ii G and let 70,71 G x M*^) be two optimal plans 

that reach the minimum in (6.2) for (i(cr,/io) and (((a,/ii), respectively. Let 7 G .^(M'^xM'^x 
M*^) be a 3-plane such that Pi, 2#7 = 70 and Pi, 3#7 = 71 where Pij denotes the projections 
on the coordinates Xi and Xj. A generalized geodesic with base point a connecting /io to 
is defined by /i^ = ((1 — t)P 2 + tP^)^^, for t G [0,1]. 

Although the theory in previous sections can be used to analyze general functionals in 
^ 2 {^^), we shall concentrate our attention in the following cases: 

The time-dependent potential energy 


V{t,fi) 


V{t, x) dy.{x), 


(6.4) 


where V : [0, 00 ) x ^ M is a time-dependent potential and the time-dependent inter¬ 
action energy 


>V(f,/i) = ^/ W{t,x,y) d{fi X fi){x,y), (6.5) 

^ JM.'ixRd, 

where W : [0, -|-cxo) x x —)■ M is an interaction potential. We also are interested in the 

case of time-dependent diffusion coefficient in the internal energy functional 


U{t,fi) = K,{t)U{y.) = K,{t) 


plog(p) dx, 


( 6 . 6 ) 


where k : [0, cxo) — )■ (0, 00 ) and dfx = p dx is an absolutely continuous measure with respect 
to the Lebesgue one. For singular measures, we set W(f, p) = -|-cxo. 
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Remark 6.2. The tools developed in previous seetions are not direetly applieable for time- 
dependent K. The reason is that the condition E3 is not satisfied for arbitrary n, but only for 
K constant. So, we postpone the case of n depending on t for later. 


6.1 The case with constant diffusion 


We consider the functionals 


Siit.n) = i4A{ii) + V(f,/i) 


and 




where U is defined in (6.6) and k > 0 is a constant. In order to apply the theory, we assume 
some conditions on the potentials. 

VI.- For each fixed f > 0, l/(f, •) is A(t)— convex, for some function A : [0, cxo) — )■ M in 
oo)), that is, V (f, x) — is convex. 

V2.- Let d°V{t, x) denote the element of minimal norm in the subdifferential of lA(f, •) at 
the point x G M'^. We assume that \d°V (f, 0)| is locally bounded and t ^ V{t, 0) is 
locally bounded from below. 

V3.- There exists a function f3 G ^L([0,+oo)) such that 



(6.7) 


We consider V2 for x = 0 only for simplicity. Indeed, this condition can be assumed for 
any (fixed) xq G M'^. Moreover, it is not necessary to choose the element of minimal norm 


in the subdifferential. In fact, it would be enough to make a measurable choice (in t) in the 
subdifferential. 


We start with the following result. 


Proposition 6.3. Assume the hypotheses VI to V3. If there exists p G such that 

dp = pdx, 



dx < oo, and / 1/(0, x) dp{x) < oo, (6.8) 


then the functional Si satisfies El to E5. 


Proof Taking s = 0 in V3, it follows that 



(6.9) 
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Next, we show the estimate 


V(t,x) > —A(t) — (6.10) 

where A{t) = —V{t, 0) + ^\d°V{t, 0)p and = |(1 + A^) for all t G [0,T]. In fact, by 
the definition of subdifferential, we have 


V{t,x) > l/(t, 0) + 0), x) — 

> v'((,o)-i|a“V'((,o)p-i(i + A^)|xp, 

and so (6.10) follows. This estimate implies that the functional V is lower semicontinuous 
with respect to the Wasserstein metric, for each fixed f > 0. Since the internal energy func¬ 
tional U is also lower semicontinuous (see [26]), we obtain El. Using (6.9), the second con¬ 
dition in (6.8), and (6.10), it follows that Dom()2(f, •)) is nonempty and time-independent, 
which gives E2. The property E3 is a direct consequence of V3 by taking = (5 o £ 

E5 follows by using the convexity of the function /r —d 2 (cr, /i) along generalized geodesics 
with base point a and the convexity of the potential V{t, •). Next, we turn to E4. Recall the 
estimate [11] 

[ plog(p) dx >-U(l-f M 2 (/i))“, (6.11) 

jRd- 

where a E (0,1) and U > 0 are constants depending only on the dimension d, and d/r = 
pdx E Thus, we obtain from (6.10) that 

= - -®) - ''C(l + M2(.p)r - A{t). 


Choosing r*(T) > 0 such that > 1 -f Aj., and using V2, the last expression is bounded 
from below by a constant depending on a, k, d, A^, r*, T, and so E4 follows. □ 


In view of the hypotheses in Theorem 5.4, we need to impose one more condition on V 
in order to obtain the needed regularity for the functional V, as expected. 


Lemma 6.4. Let k > 0 and Di = Dom{Si). If, in addition to VI, V2 and V3, we assume that 
t ^ V{t, x) is differentiable for each x E then the function t —)■ V{t, p) is differentiable 
for each p E Di. Moreover, for each sequence tn ^ t and d 2 (/in, d) —>■ 0, we have that 


fim - 

n—AcxD tn — t 


A 


V{t, x) dp{x). 


( 6 . 12 ) 


Proof We take a E <^^ 2 ,ac(®'^) and the maps and that realize the optimal transports 
from cr to pn and from a to p, respectively. Then, 


)2(tn,/^n) dn) 

tn t 


U(tn,t^) -U(t,t^-) 
tn t 


da{x). 


(6.13) 
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By [26, pag. 71], we have that t[^"(a:) —)■ t^{x) a.e. in with respect to a. Using a version 
of the dominated convergence theorem, we can take the limit in (6.13), as n —)■ oo, and 
obtain (6.12). □ 

The metric space and functionals addressed here present more structure than 

those in previous sections, where an abstract theory has been developed. So, it is natural to 
wonder if gradient flow solutions as in Definition 2.3 is related to other senses of solutions 
in,^ 2 (R'^)- In this direction, we show that the solution u associated to the functional £i is 
in fact a distributional solution for the Fo kk er-Planck equation. In the next result, we state 
precisely this fact and give some properties for u. 

Theorem 6.5. Consider the functional £i with k > 0 and potential V satisfying the assump¬ 
tions VI to V3 and the differentiability condition in Lemma 6.4. Then, given po G 
the curve p : [0, oo) —)■ t^ 2 (M'^) given in Theorem 4.4 is a distributional solution for the 
Fokker-Planck equation 

dtp = ndsp+ y ■ iyv{t,x)p), (6.14) 

with limi_j.o+ pit) = po weakly as measure. If n > 0, such curve is absolutely continuous 
with respect to the Lebesgue measure, i.e. dpt{x) = p{t,x)dx, and pit, ■) G IF;q^^(M'^). Also, 
p satisfies the energy identity 


£i{s,p{s)) = £i{t,p{t)) + 
for s < t, where \l/i : [0, oo) x —)■ 



(|T'i(r,f)p - dtV{r,x)) dpfix) dr 


is a vector field satisfying the identity 


(6.15) 


p{t, x)^'i(f, x) = nVp{t, x) + p{t, x)'VxV{t, x), for k > 0, 


(6.16) 


and = d°Vit,x) for n = 0. Moreover, if the function A satisfies E6, and pi, p 2 are two 
solutions, we have the contraction property 


d2{pi{t), P2{f)) <6-^0 ''*d2(/io, Pi)- 

Proof. First we calculate the variation of £i{t, p{t)). We have that 


(6.17) 


£i{s,p{s)) - £i{t,p{t)) = nil{{p{s))-l{{p{t))) + (U(s,a;) - U(f,a;)) 46,18) 


+ / V {t, x) dps 


V if, x) dpt- 


(6.19) 


Dividing (6.19) by s — t, using Lemma 6.4, and recalling that the function £i{t, p{t)) is 
absolutely continuous, we get 


A 

dt 


£i{t,p{t)) = 


(^i(t, x),v{t, x)) dpt{x) + / dtV{t, x) dpt{x). 


( 6 . 20 ) 
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where v : [0, cxo) x —)■ is the veetor field assoeiated to the absolutely eontinuous 

eurve Ht, •)llL 2 (^,;Kd) = and is the veetor field satisfying ||^i(t)||i 2 (^^.jRd) = 

Moreover, v verifies the eontinuity equation 

dt^it + V • = 0, (6.21) 

in the distributional sense, with Vt{x) = f (t, x). Using (5.7) together with (6.20), we obtain 
—'l/i(t, x) = n(t, x) for /ij-a.e. x G and the identity (6.15). □ 

Similar results hold true for the funetional 82 - In what follows, we state the hypotheses 
for W and W in order to treat 82 in light of the abstraet Theorem 5.4 in metric spaces. 

Wl.- For each fixed f > 0, the interaction potential W{t^x^y) is symmetric and, for 
f = 0, it satisfies a quadratic growth condition, namely W (t, x,y) = W (f, y, x) and 
W{Q,x,y) < C{1 + |xp + ||/|2). 

W2.- For each fixed f > 0, IU(f, •) is A(f)—convex, for some function A : [0, cx)) —M as 
in E5. Let d°W{t, x, y) denote the element of minimal norm in the subdifferential of 
W {t, •) at the point (x, y) R'^. We assume that \d°W (f, 0,0) | is locally bounded 

and t —)■ hF(f, 0, 0) is locally bounded from below. 

W3.- There exists a function f5 G ^L([0,+oo)) such that 

|IU(s,x, I/) — VF(t,x, |/)| < [ /3(r) dr(l + |xp + ll/p), for 0 < s < f. (6.22) 


The reason for assuming Wl is to obtain a quadratic growth for W(t, x, y) , for each f > 0, 
and then one can use the results in [5]. In fact, using Wl, this growth follows directly from 
W3. Proceeding as in Proposition 6.3, again we get that the functional 82 satisfies El to E5. 
Assuming a differentiability property in the f-variable, we obtain the analogous of Lemma 
6.4. Here we only state the results for the functional 82 . The proof is similar to that of 
Theorem 6.5 and is left to the reader. 

Theorem 6.6. Consider the functional 82 with the interaction potential W satisfying Wl to 
W3. Suppose also that for {x,y) G R'^ x R'^ the function t —?■ W{t,x,y) is differentiable. 
Then, given /tq G the curve p : [0, 00 ) —)■ ^ 2 (M'^) given in Theorem 4.4 is a 

distributional solution for the continuity equation 

dtp = S7 ■ {y{t,x)p), (6.23) 

with limt_j.o+ pit) = Pq weakly as measure, where 

v(f, x)= / x,y)p{y) dy, p-a.e. in R'^, (6.24) 
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and rjit, x, y) = x, y) + r} 2 {t, y, x)) for some Borel measurable selection (r^i, 772 ) ^ 

dW (t, •, •)• Moreover, fi satisfies the energy identity 

82 ( 3 , fi{s)) = 82 ( 1 , + [ [ {\v{r,x)\‘^ - dtW{r,x)) dfir{x) dr (6.25) 

Js Jr<^ 

for s < t. Furthermore, if the function A satisfies E6 and fii, /i 2 are two solutions, we have 
the contraction property 

d 2 {pi{t),p 2 {t)) < ‘^M2(/io,/ii)- (6.26) 

Remark 6.7. Let us observe that the vector field v in (6.23) is characterized by the form 

(6.24) thanks to the results of Carrillo-Lisini-Mainini [5]. They showed that, in general, 
the Borel measurable selection of p depends on the probability p G t^ 2 {'^'^) and is not 
necessarily given by the minimal selection in the subdifferential of W. In the particular case 
when W (t, x, y) = w{t, y — x) is given by a symmetric function w : [0, 00 ) x —)■ M, the 
A(f)-convexity of kL(t, •, •) follows from the one of w only if X{t) <0 and therefore we can 
not use the results for any A (t)-convexity of w. 

Of course, we can consider a functional of the type (see subsection below to the time- 
dependent viscosity) 

8 {t,p) =U{p) + W{t,p) 

and apply the metric theory in order to obtain existence of curves satisfying the conclusions 
in Theorem 5.4, the contraction property and a continuity equation. On the other hand, we 
do not know how to describe the velocity field v in this general case, however it is expect 
that p satisfies a Mackean-Vlasov equation of the type dtp = Ap + V ■ (v/i) for v as in 

(6.24) . 

Finally, if we assume that W{t,x,y) = w{t,y — x) is A(f)-convex with A(t) < 0 and 
satisfies a doubling condition property w{t,x + y) < Cfil + w{x) + w{y)), then one can 
show that the curve pt given in Theorem 4.4 is a distributional solution of the Mackean- 
Vlasov equation 

dtpt = Apt -f V ■ ((VM;(f) * pt)pt). 


6.2 The case with time-dependent diffusion 

Now we consider the case when k : [0, 00 ) ^ (0, 00 ). For the sake of simplicity, we 
consider the functional 

8 (t,p) = K(t)U(p) + V(t,p), (6.27) 

where Id and V are defined in (6.6) and (6.4), respectively, and k is a positive function locally 
absolutely continuous. Also, we assume that VI to V3 hold true. Thus, by assuming that 
1/(0, •) satisfies (6.8), we have that the domain of 8 is time-independent. Notice that the 
functional 8 satisfies El, E2, E4 and E5, but not E3. In fact, as observed in Remark 6.2, E3 
holds true, if and only if, K(t) is constant. Here we need to assume that k is non-increasing. 
An important fact is the following: 
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Remark 6.8. Let fin £ and G [0, cxo) be two bounded sequences, where fin is 

bounded with respect to the Wasserstein metric d 2 , such that the numeric sequence 8/i„) 
is bounded from above. Then, the numeric sequences U{fin) and A^n) are bounded. In 
fact, it follows from (6.11) that the sequence K{tn)U{fin) is bounded from below and thus 
V(fn, /^n) is bounded from above. Similarly, from (6.10) we have that V(fn, k^n) is bounded 
from below, and then K,{tn)U{jjin) is bounded from above. 

Using Remark 6.8, we obtain easily the same conclusions of Lemmas 3.1 and 3.2. Let the 
potential V be differentiable in the f-variable. Using the minimality of we have that 
for To < Ti 

ri(u) - (ft+ro,To(u) < {K{t + Ti) - K{t + Tq))W(/ it^"'°) + V{t + Ti, <+"'“) 

- V{t + To, 4+"°) 

< {K{t + Ti) - K{t + ro))W(Aitr“) + 

rt+Ti 

+ / P{r) dr{l + M{ul+^°)), 

J i+ro 

where above we used the estimate (6.11). Analogously, the reverse inequality follows. In 
view of Remark 6.8, we can argue as in Proposition 3.4 and Corollary 3.5 in order to obtain 
the identity (3.8) for the functional (6.27). 

Up until this point, notice that we have not needed the monotonicity hypothesis for n. 
In what follows, we comment on an essential step in order to recover Lemma 3.8. In fact, 
recalling the notation for discrete solution (Ul) of the variational scheme, and using that K{t) 
is non-increasing and (6.11), we estimate 

n 

Y^4timui-')-u{ui)) < ^mwi) - 

n 

-cY^wi.) - K(«i“‘))(i + w(£/r‘)). 

i=i 

From here we can repeat the arguments in order to obtain the same conclusion of Lemma 
3.8. In the case when 12 = 0, it is not necessary to suppose the monotonicity of n because 
the difference W(f/^“^) —U{Ul) is positive and a more direct estimate can be performed. 
Going back to Section 4, it is easy to see that it remains only to estimate 

[(1 - - K{rAs)))U{U^{s))]+ ds, for 0 < f < T, (6.28) 

where T > 0 is fixed and r, rj are two partitions with small sizes. Indeed, since 
S^T-rit), U-rit)) is bounded from above by a constant independent of r, it follows from Re¬ 
mark 6.8 is bounded by a constant independent of r. Thus, the integral (6.28) 
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can be estimated by proceeding similarly to Proposition 4.3, and then we obtain the con¬ 
vergence of the approximate solutions (4.8). In this way, the functional £ defined in (6.27) 
presents properties and results contained in Sections 4 and 5. So, we have the following: 

Theorem 6.9. Let £ be the functional defined in (6.27) with k : [0, cxd) —)■ (0, oo) an 
non-increasing absolutely continuous function and let the potential V satisfy VI to V3 
and the differentiability condition in Lemma 6.4. Then, given pq G the curve 

p : [0, oo) —)■ obtained in Theorem 4.4 is absolutely continuous with respect to the 

Lebesgue measure, i.e. dpfix) = p{t,x)dx, p{t, •) G (Ef) for each t G [0, cxo), and p is 
a distributional solution for the Fokker-Planck equation 


dtp = K,{t)Ap -f V ■ (fi7V{t,x)p), 

with limt_j.o+ pit) = po weakly as measure. Also, pft) satisfies the energy identity 


(6.29) 


£iis,pis)) = £iit,pit)) -f 



(|4/i(r,f)|^ — dtVir, x))pir, x) dx dr 



where 4/i : [0, oo) x 


-)■ 


nfr)pir, x) log(p(r, x)) dx dr, for s < t, 


is a vector field satisfying 


pft, a;)4^i(f, x) = K(f)Vp(f, x) + pft, x)W{t, x) for pt-a.e. x G 


(6.30) 


(6.31) 


Moreover, if the function A satisfies E6 and pi, p 2 are two solutions, we have the contraction 
property 

d2(/ii(f),/i2(f)) < '^M2(/io,/ii)- (6.32) 


6.3 More general internal energy 

In this subsection we give the outline to construct the time-dependent gradient flow for 
more general internal energy functionals. Let U : [0, oo) x [0, oo) —?■ M be a continuous 
function such that C^((0, oo) x (0, oo)). Consider the internal energy functional 


j !m.<iUit,pix)) dx, if dp = pdx 
\ -foo, otherwise. 


(6.33) 


We assume the following condition on U. 


Ul.- There exist functions a, A : [0, oo) ^ [0, oo) with a e L,^„^([0,oo)) and A G 
L^([0, oo)) such that 


dU 

Ait)U^f),z) <—it,z) <ait)U (0,z), 


(6.34) 


for all t,z G [0,-foo), and f/(0,z) has superlinear growth at infinite, i.e. 

li„, LM = +00, 
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U2.- There exist a E (0,1) with a > ^ and positive eonstants ci, C 2 > 0 sueh that 


f^( 0 , z) > —Ciz — C 22 ". 


U3.- U(0,0) = 0, z —>■ U (t, z) is eonvex, and z —)■ z'^U (t, z is eonvex and non-inereasing 
on ( 0 , + 00 ), for eaeh t > 0 . 

Without loss of generality, we ean assume that ||A||i = A(t) dt < 1; otherwise, we 
ean replaee U by pjj^- Firstly, let us note that U1 and U2 imply 

r* dU 

U{t,z) = / -—{r,z)dr + U{0,z) 


dt 


> 


A{r) dA U+{0,z) + U{0,z) 


= A{r)drjU+{0,z)-U-{0,z) 

> —U~{0, z) > —Ciz — C 22 ;". 


(6.35) 

(6.36) 


Then, reealling that a > it follows from (6.36) that 


dt(t, - I Cl + C 2 / p(a;)“ dx 


> 


Cl + C 2 ( / (1 + \x\^)p{x) dx 

= ~ (ci + 026 * 0,(1 + M 2 (/i))“) . 


(1 + |xP) i-“ 


dx 


1—Q ' 


(6.37) 


Therefore, the funetional in (6.33) is well-defined from [0, +cxo) x to (—cxo, +cx)]. 

It follows from (6.35) that U{t, ■) has a superlinear growth, for eaeh fixed t > 0. So, by stan¬ 
dard arguments (see [17]), one ean show that the funetional U{t, •) is lower semieontinuous 
with respeet to the weak topology, for eaeh fixed t > 0. Thus, U{t, •) verifies El. 

Let p E t^ 2 ,ac(^^) be sueh that dp = pdx and W(0, p) < 00 . We have 


U{t,p{x))<yJ a{r)drjU (0,p(x)) + f/(0,p(x)), 

and then p) < 00 for all f > 0. On the other hand, if dp = pdx is sueh that Z7(f, p) < 
-l-cxo for all t > 0, then, by substituting 2 ; = p{x) in (6.35), we get 


— j A{r) dr^ p{^)) ~ U (0, p{x)) < U (t, p{x)). 
It follows by integrating (6.38) that W(0, p) < 00 . So, U{t, p) verifies E2. 


(6.38) 
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Denote Dom{U{t, •)) = D C ^ 2 ,ac{^'^)- Sinee U (0, 0) = U~^{0, 0) = 0, note that 



and then D is nonempty. For s < t and /r G D with d/r = pdx, we have 



The same arguments used in (6.37) lead us to 



U{t, /i) — U{s, /i) < 


a(r) dr j (ci + + M 2 (/i))), for all 0 < s < f. (6.39) 


We are going to use (6.39) as a substitute for the eondition E3. Also, E4 follows from (6.39). 
In fact, 



Now, it is easy to see that for r* > 0 small enough the last expression is bounded from below, 
as desired. Note that we have E5 with A = 0 because U3 implies that W(f, •) is convex along 
of generalized geodesics. 

Let us remark that Lemmas 3.1 and 3.2 can be proved by proceeding as in Section 3 (and 
using (6.35) for Lemma 3.2). In order to recover the differentiability property in Proposition 
3.4, we recall the notation 


2 r 




+ W(f, z/) 




Then, by taking tq < n and /r G D, we have 



(6.40) 


Denoting = p^^'^dx, we can estimate 



The last integral on is uniformly bounded in tq on compact sets of (0, r*]. Replacing the 
roles of To and ri, we get 


— <^t+TO,To{f^) >ld{t + Tl, —U{t + TQ, 



(6.41) 
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and 


U{t + ri, - U{t + To, /xtr) > - A{r) dr [ [/+(0, pt^) dx. 

Now we need a uniform estimate for U~^{0, {x))dx. Substituting t = t + ti and 

z = (x) in (6.35), and afterwards integrating it, we arrive at 

{1-\\A\U) [ U^{0,pl^^^{x))dx< [ U-{0,pl^^^{x))dx + <^t+r,,rAf-^). (6.42) 

jRd. jRd 

The first term in the right hand side of (6.42) is loeally uniformly bounded in (0, r*]. By the 
eontinuity of the map r <o)+T- i-(p), the seeond term also verifies so. Therefore, we eon- 
clude that the function r —)■ St+r^rip) is absolutely continuous in each compact subinterval 
of (0, T*]. Now a version of the dominated convergence theorem leads us to the formula 


for each differentiability point r G (0, r*]. The identity above implies the integral equality 
(3.8) in Corollary 3.5. 

In the sequel, we sketch the proof of Lemma 3.8 in the case of this present section. Re¬ 
calling the notation for the discrete solution in (2.8), we have 


;(M,((7;) - AMU;)) < J2 - M.iUp^)) 




< 


i=i 


2Tj 




-(M(o, f/J) - XM,(ui 


i=i 




i=i 


Using the above estimate and (6.39), we can proceed as in Lemma 3.8 and reobtain the 
conclusions of this lemma for the functional W(f, p). 

Now we deal with the convergence of the approximate solutions. In comparison with 
subsection 4.2, there is only a new term that reads as 

[\i-iAt))m%{t),uAt))-i^{%{t),uAm dt. 

Jo 

Notice that it is necessary to consider only the case T-rit) < 7^(t). So, we have that 

U{TAt),UAt))-ld{%{t),U^{t))<ir"^\{r)dr] [ U^{0,U^{t,x))dx. (6.44) 

yr-rit) J Jm^ 
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Using (6.35) we can estimate the integral over in (6.44) locally uniformly in [0,cxo). 
Now, by replacing the function /3 by a or A, one can repeat the same arguments in the 
proof of Proposition 4.3, obtain the estimate (4.19) and afterwards the convergence of the 
approximate solutions (4.8) to a curve /i : [0, oo) —)■ 

In what follows, we sketch the main arguments for the convergence of the De Giorgi 
interpolation (see Proposition 5.2). 

Let 5 = t — for t G f”]. Using the minimizer property of Ut and Ut, we can 
obtain 

+ M(i. UAt)) - mK, UAt)) 

< ^^d2(U,(f),U,(f)) +W(t,I7,(t)) -W(C,U,(f)), (6.45) 


and so 


+ / a{r) dr / U~{0,UT{t,x)) dx 


+ / A{r) dr / U~^{0,Ur(t,x)) dx. (6.46) 

Jt jRd 


By using these estimates, one can obtain the conclusions of Proposition 5.2. In order to 
recover some properties in Theorem 5.4, note that the inequality > in (5.7) follows from the 
same arguments and Patou’s Lemma. The absolutely continuity of the map f —> W(f, /i(f)) 
in each bounded interval of [0, oo) is a consequence of (5.4) with \{t) = 0 and the estimate 
(6.35). 

Although we are not able to obtain the reverse inequality < in (5.7), and so the energy 
identity (5.7) (see Remark 6.12 below), it is not hard to show the estimate 

W(f,/i(f))-W(s,/i(s)) < ^ \fi'\^{r) dr 

|c)Z7(r)|^(yu(r)) dr. 

Even without the energy identity and minimal selection, it is possible to show that the 
approximate solutions converge to a distributional solution /i. Denote P(t, z) = z^{t, z) — 
U{t, z). We summarize the results for the functional (6.33) in the theorem below. 



Theorem 6.10. Consider the internal energy functional defined in (6.33) and assume U1 to 
U3. Then the curve /i : [0, cxo) d^it = p{x,t)dx, given by Theorem 4.4 is a 

distributional solution of the equation 


dtp - Vx • pit, x))) = 0 


(6.47) 


with initial condition pq G Moreover, the contraction property (5.17) holds. 
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Proof. For simplicity take a uniform step size r > 0, eonsider the partition {0 < r < 2r < 
3r < • • ■ }, and ehoose ^ G R'^). Consider also the flow assoeiated to the field 

i.e. 




(6.48) 


Then, by the minimizer property of Uf, we have 

di{ur\us) 


U{nT, Us) 


-U{nT,U:)- 


dl{Ur\Uf) 


T 


> 0 , 


(6.49) 


where Us = ^si^U^ is the push-forward of Uf via $ 5 . Then, by standard arguments (see 
e.g. [ 11 ]), it follows that 


U{nT,Us)-U{nT,Uf) 

hm-;- 

(S^-o S 


—P{nT, Uf{x))di\^ dx 


and 


ii,„ --1 Us) - i/;‘) 

( 5->-0 5 


- — —■ii.y) d-i{x,y), 


where 7 G r(f/" (7”) is an optimal plane for the transport from Uf ^ ioUf. Changing ^ 

by in (6.48) (by symmetry in (6.49)) and taking ^ = VC, we obtain that 


{x-y) 


T 


■ VC( 2 /) d-f{x,y) 


P(nr, Uf{x))AC dx = 0. (6.50) 


Let us remark that the above ealeulations also allow to eonelude that P{nT, Uf) G kF^’^(R'^) 
is bounded uniformly. Thus, we ean use an argument of weak eonvergenee and estimates as 
in Lemma 3.8 in order to obtain that the eurve y : [0, cxo) —)■ ^20^^) solves (6.47) in the 
distributional sense. □ 


Remark 6.11. The eonditions U1 to U3 work well if we eonsider a funetional as being the 
sum of the internal energy and another funetional as in two previous subseetions. In the 
present subseetion, we have preferred to eonsider only the internal energy for the sake of 
simplieity. 

Remark 6.12. Let us observe that the energy identity was not obtained in Theorem 6.10. The 
reason is that, in order to obtain such property in this general case, it would be neeessary to 
handle the limit 

lim 

t — tg 

By making a ehange of variable (see [26, Theorem 4.8]), the ealeulus of (6.51) is related to 
stability results for the Monge-Ampere equation. However, as far as we know, sueh results 
are available in the literature (see [9]) under restrietions stronger than the ones that we have 
in our eontext. 


38 









References 


[1] Agueh, M., Existence of solutions to degenerate parabolic equations via the Monge- 
Kantorovich theory, Adv. Differential Equations 10 (3) (2005), 309-360. 

[2] Ambrosio, L., Gigli, N., Savare, G., Gradient flows: in metric spaces and in the space 
of probability measures, Birkhauser, 2005. 

[3] Ambrosio, L., Minimizing movements. Rend. Aeead. Naz. Sei. XL Mem. Mat. Appl. 
(5) 19 (1995), 191-246. 

[4] Carrillo, J. Ferreira, L. Preeioso, J., A mass-transportation approach to a one dimen¬ 
sional fluid mechanics model with nonlocal velocity, Adv. Math. 231 (1) (2012), 306- 
327. 

[5] Carrillo, J. A. Lisini, S. Mainini, E., Gradient flows for non-smooth interaction poten¬ 
tials, Nonlinear Analysis 100 (2014), 122-147. 

[6] Carrillo, J. A., MeCann, R. J., Villani, C., Contractions in the 2-Wasserstein length 
space and thermalization of granular media, Areh. Ration. Meeh. Anal. 179 (2) (2006), 
217-263. 

[7] De Giorgi, E., New problems on minimizing movements. Boundary Value Problems for 
PDE and Applieations, C. Baioeehi and J. L. Lions, eds., Masson, 1993, 81-98. 

[8] De Giorgi, E., Marino, A., Tosques, M., Problems of evolution in metric spaces and 
maximal decreasing curve, Atti Aeead. Naz. Lineei Rend. Cl. Sei. Fis. Mat. Natur. (8) 
68 (1980), 180-187. 

[9] De Philippis, G., Figalli, A., Second order stability for the Monge-Ampere equation 
and strong Sobolev convergence of optimal transport maps. Anal. PDE 6 (4) (2013), 
993-1000. 

[10] Evans, L. C., Gariepy, R. E, Measure Theory and fine properties of functions. Studies 
in Advanced Mathematics. CRC Press, (1992). 

[11] Jordan, R., Kinderlehrer, D., Otto, E, The variational formulation of the Fokker-Plank 
Equation, SIAM J. Math Anal. 29 (1) (1998), 1-717. 

[ 12] Jun, Chanyoung, Pursuit-evasion and time-dependent gradient flow in singular spaces. 
Thesis (Ph.D.)-University of Illinois at Urbana-Champaign, 2012. 76 pp 

[13] Lisini, S., Nonlinear diffusion equations with variable coefficients as gradient flows in 
Wasserstein spaces, ESAIM Control Optim. Calc. Var. 15 (3) (2009), 712-740. 

[14] A. Lytchak, Open map theorem for metric spaces, St. Petersburg Math. J., 16, (2005), 
1017-1041. 


39 


[15] Marino, A., Saccon, C., Tosques, M., Curves of maximal slope and parabolic varia¬ 
tional inequalities on nonconvex constraints. Ann. Scuola Norm. Sup. Pisa Cl. Sci. 16 
(2) (1989), 281-330. 

[16] U. F. Mayer, Gradient flows on nonpositively eurved metrie spaces and harmonic maps. 
Comm. Anal. Geom., 6, (1998), 199-253. 

[17] McCann, R. J., A convexity principle for interacting gases, Adv. Math. 128 (1) (1997), 
153-179. 

[18] Mielke, A., Rossi, R., Savare, G., Nonsmooth analysis of doubly nonlinear evolution 
equations. Calc. Var. Partial Differential Equations 46 (1) (2013), 1-2, 253-310. 

[19] Mielke, A., Rossi, R., Savare, G., Variational convergence of gradient flows and rate- 
independent evolutions in metric spaces, Milan J. Math. 80 (2) (2012), 381-410. 

[20] Otto, R, The geometry of dissipative evolution equations: the porous medium equation, 
Comm. Partial Differential Equations 26 (1-2) (2001), 101-174. 

[21] Petrelli, E., Tudorascu A., Variational Principle for General Diffusion Problems, Appl. 
Math. Optim. 50 (3) (2004), 229-257. 

[22] Roche, T., Rossi, R., Stefanelli, U., Stability results for doubly nonlinear differential 
inclusions by variational convergence, SIAM J. Control Optim. 52 (2) (2014), 1071- 
1107. 

[23] Rossi, R., Mielke, A., Savare, G. A metric approach to a class of doubly nonlinear 
evolution equations and applications, Ann. Sc. Norm. Super. Pisa Cl. Sci. 7 (1) (2008), 
97-169. 

[24] Veretennikov, A. Yu, On ergodic measures for McKean-Vlasov stochastic equation, 
Monte Carlo and quasi-Monte Carlo methods 2004, 471-486, Springer, Berlin, 2006. 

[25] Villani, C., Optimal transport: old and new, Grundlehren der Mathematischen Wis- 
senschaften [Eundamental Principles of Mathematical Sciences] 338, Springer-Verlag, 
Berlin, 2009. 

[26] Villani, C., Topics in Optimal Transportation, Graduate Studies in Mathematics 58, 
American Mathematical Society, providence, RI, 2003. 


40 


